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Creatures as varied as mammals, fish, insects, reptiles, and birds have an intrigu- 
ing 'sixth' sense that allows them to orient themselves in the Earth's magnetic field. 
Despite decades of study, the physical basis of this magnetic sense remains elu- 
sive. A likely mechanism is furnished by magnetically sensitive radical pair reactions 
occurring in the retina, the light-sensitive part of animal eyes. A photoreceptor, 
cryptochrome, has been suggested to endow birds with magnetoreceptive abilities 
as the protein has been shown to exhibit the biophysical properties required for an 
animal magnetoreceptor to operate properly. Here, we propose a theoretical analy- 
sis method for identifying cryptochrome's signaling reactions involving comparison 
of measured and calculated reaction kinetics in cryptochrome. Application of the 
method yields an exemplary light-driven reaction cycle, supported through transient 
absorption and electron-spin-resonance observations together with known facts on 
avian magnetoreception. The reaction cycle permits one to predict magnetic field 
effects on cryptochrome activation and deactivation. The suggested analysis method 
gives insight into structural and dynamic design features required for optimal de- 
tection of the geomagnetic field by cryptochrome and suggests further experimental 
and theoretical studies. 



E-mail: ilia@illinois.cdu 
E-mail: kschulte@ks.uiuc.edu 



2 



I. INTRODUCTION 



Migratory birds travel annually thousands of kilometers, navigating by using various cues, 
including the Earth's magnetic field jl-5|; non-migrant bird species utilize the geomagnetic 
field similarly to find their way back to the breeding nest |6l-|l0| . But how do birds de- 
tect the geomagnetic field? This question focuses on one of the most fascinating unsolved 
mysteries of sensory biology. At a first glance, tiny iron-oxide particles detected in the up- 
per beak of some bird species llMlq provide a natural explanation for the avian magnetic 
compass sense, but behavioral studies revealed that the ability of night-migrating birds to 
perform magnetic compass orientation is affected by the ambient light 17H2lj, which does 



not penetrate through the beak skin. The latter observation lead to the suggestion that a 



photochemical reaction in the bird's eyes produces spin-correlated radical 



the sensor embodying a magnetic field- dependent signaling cascade ! 



22 



pairs that act as 



27| . This sugges- 



tion is based on the observation that the recombination reactions of spin-correlated radical 
pairs, can be magnetic field-sensitive 



28 



29| . The radical pair hypothesis gained stronger 



support when it was experimentally demonstrated that the magnetic compass in birds is 
still functioning properly after the iron-mineral-based receptors in the beak are deactivated 
3o|-33|. 

Spin-correlated radical pairs are typically created in a photochemical reaction that pro- 
duces radicals from a molecular precursor in either an electronic singlet (S) or triplet (T) 
state. Under the infiuence of intramolecular electron-nuclear hyperfine interactions, the rad- 
ical pair oscillates between the S and T states, a process known as S O T interconversion. 



External static magnetic fields affect the rate of this process and, hence, a. 
the respective reaction products formed from the S and T radical pairs 29 



ter the yields of 
. During 



34 



the last years the radical-pair hypothesis of the avian magnetic compass gained significant 



experimental 40l-l44j and theoretical 22 



-2) 



45 



46[ support. The first experimental proof- 



of-principle, demonstrating that under a static magnetic field, as weak as that of the Earth, 
a chemical reaction can act as a magnetic compass by producing detectable changes in the 
chemical product yield, was achieved for a carotenoid-porphyrin-fuUerene model system 47 |. 

However, for the radical-pair mechanism to play a role in magnetoreception, molecules 
with certain biophysical characteristics must exist in the eyes of migratory birds. Cryp- 



tochromes 



48 



-|51|. a class of photoreceptor proteins, were proposed as the host molecules 
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for the crucial radical pair co-factors that putatively act as a primary magnetoreceptor 23 1. 

ml 



Cryptochrome is a signaling protein found in a wide variety of plants and animals |48 
Its role varies widely among organisms, from the entrainment of circadian rhythms in ver- 
tebrates to the regulation of hypocotyl elongation and anthocyanin production in plants 
52l-l54| ■ It has recently been demonstrated, using UV/visible transient absorption and 



electron paramagnetic resonance spectroscopy |55l. l56l|. that vertebrate cryptochromes form 
long-lived radical pairs, involving a flavin radical and a radical derived from a redox-active 
amino acid. This observation demonstrates that cryptochrome harbors the type of radical 
pair needed for magnetic compass action and suggests that cryptochromes, through S ^ T 
interconversion, are influenced by an external magnetic field. Radical pairs seen in DNA 
photolyase 57H63l|. a light sensitive DNA repair enzyme closely related to cryptochrome, 



support the observation 



55, 



56j in cryptochrome; indeed a recent study 59|] demonstrated a 



magnetic field effect in DNA photolyase, suggesting a similar effect to exist in cryptochrome. 

It has been verified that cryptochromes exist in the eyes of migratory birds 55|, l64l-l67l|. 
and that at least some cryptochrome-containing cells within the retina are active at night 



when the birds perform magnetic orientation in the laboratory |64j . I67l |. Furthermore, a 
distinct part of the forebrain, which primarily processes input from the eyes, is highly active 
at night in ni ght-mig ratory garden warblers {Sylvia borin) and European robins {Erithacus 
rubecula) 33|, IGSMTII]. In summary, many findings are consistent with the hypot hesis that 



magnetic compass detection in migratory birds takes place in the eye [33|, l72N76| and that 
cryptochromes are the primary magnetoreceptors. However, despite the success of these 
findings, a completely satisfactory description of the mechanism of the magnetic field effect 
in cryptochrome is still missing. 

The process of cryptochrome photoactivation was discussed earlier and several reaction 



schemes were proposed 



57, 



66 



77|-8^. However, the proposed schemes are incomplete, usu- 



ally accounting for only some of several observations. Cr ypto chrome binds internally the 



chromophore flavin adenine dinucleotide 



fFAD) 



tochromes from Amhidopsis thaliana 55 



77-79 



50 



66 



77 



79 



-|82|. At least in plant cryp- 



84| blue light leads to conversion of the 



fully oxidized FAD to the semireduced FADH* form, the latter representing the signaling 
state. The conversion happens in the course of a light-induced electron transfer reaction 
involving FAD and a chain of three tryptophan amino acids that bridge the space between 



FAD and the protein surface 



55 



56 



85|, as illustrated in Fig. [TJ Atomic level structures of 
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plant cryptochrome are known for Arabidopsis thaliana cryptochrome-1 8l| and Arabidopsis 



thaliana DASH-type cryptochrome-3 



82l |. In cryptochromes from insects, light excitation 



leads to formation of a flavin anion radical, FAD" 



84 



It is currently under debate 



whether the anion radical represents the signaling state or whether it is only a function- 



ally insignificant short-lived intermediate 80|, l86|, |87j. The difference in the photocycles of 



plant and insect cryptochromes illustrates that detailed analysis of the transient states in 
cryptochromes is needed in light of the fact that major variations arise between different 
cryptochromes, even regarding the nature of their signaling states. Such analysis can only 
be performed in vitro on cryptochromes extracted from specific organisms. 

Unfortunately, relative to plant and insect cryptochromes, little is known at present 
about structure and photocycle of avian cryptochromes. Therefore, we study here a model 
"cryptochrome" from a general perspective. Our goal is to demonstrate what observations 
are needed for resolving photoactivation and magnetoreception of the protein. By suggest- 
ing an analysis method that deduces from experimental observations a reaction scheme, it 
is our hope that further physico-chemical experiments will be initiated on cryptochrome 
photoactivation that ultimately pinpoint the complete photoreaction and signaling process 
in cryptochrome and, in particular, establish when and why the anionic, FAD*~, and the 
neutral, FADH*, radicals are formed. 

One can safely state that, in spite of numerous experimental observations of intermediate 
states in cryptochrome, the actual photoreaction is still not fully understood 



77 



79 



55 



H57I, 



66 



91| . Figure |2] shows a scheme of cryptochrome activation and inactivation. 



The scheme incorporates the key observation that the flavin cofactor in cryptochrome is 



observed in three interconvertible redox forms, FAD, FADH', and FADH' 



89 



56, 3, 



77 



83 



92|. In this scheme the FAD form is inactive (non-signaling) and accumulates to a 
high level in the dark 77H80L l88| . Blue light triggers photoreduction of FAD to establish 
a photoequilibrium that favors FADH* over FAD or FADH~ 77H80l. |88| . As pointed out 
above, in plant cryptochromes, signaling is linked to formation of the FADH* state. This 
state can absorb a second, blue-green light photon, in which case FADH* is converted to the 



fully reduced, inactive 



resting state [25 



57 



59 



brm FADH : the latter reoxidizes in the dark to the original FAD 



77 



89 



9j. 



The tryptophan triad in cryptochrome, depicted in Fig. [H is crucial for the functioning 



of the protein as a primary electron donor 



55 



56 



78| . In Arabidopsis thaliana the triad 
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Trp377 
Trp324 Trp400 




FIG. 1. Arabidopsis thaliana cryptochrome in transient absorption experiment. (Top) 



8l[ is shown together with the high- 



The structure of cryptochrome- 1 from Arabidopsis thaliana 
lighted flavin cofactor (FAD) and the tryptophan triad Trp400, Trp377, Trp324. (Bottom) Tran- 
sient absorption of cryptochrome is probed by means of the pump-probe experiment, as described 
in experiment 



551]. The sample containing cryptochrome is irradiated with a pulsed laser beam 
(red) that generates a measurable concentration of excited states (FAD*) in the system. FAD* 
decays then in a series of intermediates back to FAD, while some of the intermediates being probed 
by the probe beam (white) used to measure the absorption spectrum of transient species. Ad- 
ditionally, magnetic field effects in cryptochrome can be studied if the sample is subjected to an 
external magnetic field. 
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FIG. 2. Cryptochrome activation and inactivation reactions. Cryptochrome is activated 
through absorbing a blue-light photon by the flavin cofactor, responsible for protein's signaling. 
Initially, the flavin cofactor in cryptochrome is present in its fully oxidized FAD state. After 
absorbing a photon, FAD becomes promoted to an excited FAD* state. FAD* is then protonated 
and receives an electron from a nearby tryptophan (see Fig. [T]), leading to the formation of the 
[FADH* + Trp*] radical pair, which exists in singlet and the triplet overall electron spin states, 



denoted as ^ [ • • • ] and ^J_;_ • ^ respectively. The Trp* radical 
from a nearby tyrosine 



551, 



24 



may receive an additional electron 
quenching the radical pair and 



78l |. or become deprotonated 
fixing the electron on the FADH* cofactor. Under aerobic conditions, FADH* slowly reverts back 
to the initial inactive FAD state through the also inactive FADH" state of the flavin cofactor. 
Before the [FADH* + Trp*] radical pair in cryptochrome is quenched, the electron of the FADH* 
radical can back-transfer to the Trp* radical, thereby also ending the signaling state of the protein. 
The electron back-transfer leads to the formation of FADH"*" and can only occur if the spins of the 
two unpaired electrons in the radical pair [FADH* + Trp*] are in an overall singlet state. While the 
flavin cofactor is in its FADH* state, cryptochrome may absorb a second blue-green photon, thereby 
transferring it into a non-signaling state with the flavin cofactor in the FADH" conformation. 
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consists of Trp324, Trp377, and Trp400 as shown in 
warbler {Sylvia borin) |55|, Drosophila melanogaster 



ig. jH In cryptochromes from garden 



56 



79 



80 1 and Homo sapiens jSg. 



the tryptophan-triad is conserved 24. 155l-l57l|. Before hght activation, the flavin cofactor of 



cryptochrome is present in its fully oxidized FAD state (see Fig. [2]). FAD absorbs blue light, 
being thereby promoted to an excited state, FAD*, which is then protonated, likely from a 



nearby aspartic acid 9J], and receives an electron from the nearby tryptophan (see Fig. [T]), 
leading to the formation of a [FADH* + Trp400*] radical pair 2J| . This radical pair is further 



transformed into a [FADH*+Trp324*] radical pair (denoted in Fig.|2]as ^'^[FAD 



H'+Trp*]) via 



sequential electron transfer, also ca^ 
the tryptophan triad chain 24, [sgI, 



95|, involving 



led paramagnetic-diamagnetic exchange 

. The [FADH* + Trp*] radical pair exists in singlet 
and triplet states, denoted in Fig. |2]as ^[FADH* + Trp*] and '^[FADH* + Trp*], respectively. 

The [FADH* + Trp*] radical pair is associated with the cryptochrome signaling state, 
since it involves the flavin cofactor in the FADH* redox state. We note, though, that the 
semiquinone FADH* state of FAD was identified as the signaling state in plant cryptochrome 



from Arahidopsis thaliana 77], l83|, , but that the signaling state of avian cryptochrome is 



still unknown. In our study we generalize the observation in Arahidopsis thaliana to birds; if 
future studies will demonstrate that a different redox state of the flavin cofactor is governing 
the signaling behavior of avian cryptochromes, the suggested scheme can be readily adapted. 
We also note that the Trp* radical may receive an additional electron from a nearby tyrosine 



55 



78| . or become deprotonated 2^, l93|, stabilizing the FADH* redox state. The charge 



state of the tryptophan radical in the radical pair state might be crucial for the energetics 
and kinetics of the reaction pathway, but it is not essential in a methodological sense, the 
main focus of the present study. The lifetime of the radical pair state [FADH* + Trp*] in 
cryptochrome is > 6 yus, according to transient EPR measurement {sgI . 

Under aerobic conditions, the FADH* redox state reverts back to the initial FAD state 



77 



92| (see Fig^^^ 



millisecond time scale 



55 



superoxide radical 



flavin cofactor 



25 



77 



his process is not well understood, but seems to occur on a 



It has been suggested that the back-reaction involves the 



0^(43], the reaction evolving through the inactive FADH state of the 
43L l96| , as depicted in Fig. |2l The lifetime of the FADH" state, populated 



such, should be long enough to allow its detection in transient absorption measurements 



55 



56]. 



Quenching of the [FADH* + Trp*] radical pair in cryptochrome can also arise through 
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electron back-transfer from FADH* to a tryptophan. This back-transfer, leading to formation 
of FADH"*" (see Fig. [2]), can only occur if the spins of the two unpaired electrons are in an 
overall singlet state. An external magnetic field can influence the overall electron spin state 
through the Zeeman interaction acting jointly with hyperfine coupling with hydrogen and 



nitrogen atoms 24 



25 



37 



45| . If the overall electron spin state of [FADH* + Trp*] is 



triplet, electron back-transfer and formation of FADH+ cannot occur, extending the time 
cryptochrome stays in its signaling state. The FADH"*" state has a short lifetime and decays 
quickly to the fully oxidized FAD configuration via deprotonation. 



According to recent experimental measurements [57|, |77|, |79|, |88| and as pointed out above, 
cryptochrome may absorb an additional blue-green light photon, thereby transferring itself 
into a non-signaling state, namely with the flavin cofactor in the FADH~ redox state. The 
process involves the transition FADH* — )■ FADH** — )■ FADH~, depicted in Fig. El After 
excitation, the flavin cofactor in the FADH** conformation has an electron vacancy, allowing 
it to accept an electron from a nearby tryptophan and, thereby, be transformed into FADH^. 
In order to stabilize FADH", the terminating tryptophan of the triad should be in its radical 
state. In the absence of such tryptophan radical prior to FADH* excitation, the excited 
electron in FADH** will back-transfer to the electron vacancy in tryptophan as soon as 
FADH" is formed. 



If cryptochrome is indeed the avian magnetoreceptor protein, its signaling state is ex- 
pected to be sensitive to a magnetic field, the latter controlling the lifetime of the FADH* 
state. In earlier studies we demonstrated that weak magnetic fields, comparable with the ge- 



omagnetic field, can alter the lifetime of FADH* by 5%-10% |24l. |25| . In these earlier studies 
we considered reactions in isolation, while in the present study we investigate the complete 
activation cycle. In the following we suggest a light-driven reaction cycle in cryptochrome 
that exhibits a representative magnetic field effect on the signaling state of the protein. The 
reaction cycle is supported by transient absorption and electron-spin-resonance observations 



55 



56| and incorporates known attributes of avian magnetoreception. We demonstrate how 



the proposed reaction cycle can be compared to experiment. 
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II. METHODS 



Cryptochrome is a sensory protein [56|, |66|, l77|, l79|, l88|, |89[ . Its relevant properties, namely, 



light activation and deactivation, arise through a complex reaction scheme linking many in- 
termediate states. In case of magnetoreception, cryptochrome activation and deactivation is 
apparently magnetic field dependent. The long-range goal of our study is to establish cryp- 
tochrome's photo-reaction scheme unequivocally, combining measurement and theoretical 
analysis. 



Key for the stated goal are measurements of cryptochrome transient photoabsorption 55 



561] that presently suggest different reaction schemes, consistent with the same experimental 
observations. For the sake of concreteness we single out one reaction scheme, namely the 
one shown in Fig. [2j We will demonstrate how this scheme can be reconciled with observed 
time-dependent absorption spectra. The method can be applied to variant schemes, too. 

The efficiency of light absorption at wavelength A by an absorbing medium is characterized 
through the absorbance A{X) defined as jo^] 

N 

AiX)=dY,e^{X)c„ (1) 

i=l 

where is the number of different light-absorbing components in the system with the con- 
centrations Cj (expressed in mol ■ L~^) and the molar absorption (extinction) coefficients Si{X) 
(in L ■ mol~^cm~^); d is the thickness of the absorbing medium (in cm). The molar absorp- 
tion coefficients in Eq. ([T]) are wavelength-dependent, resulting in a wavelength-dependence 
of the sample absorbance. The time evolution of the intermediate states concentrations can 
be calculated from the set of coupled kinetic equations 

=-k, [FAD] + ki'J, [FAD*] + [FADH+] + k^^J [FADH"] , (2) 
d[FAD*] ^ ^^p^^^ _ _ fc(iJ[FAD*], (3) 



dt 

d[HFADH'Trp')] ^^^^^p^j^.^ _ j^^(i) + + + [i(FADH-Trp-)] + 

fc,/[3(FADH'Trp')] + A;fJ[i(FADH'*Trp')], (4) 

d[3(FADH'Trp')] ^ ^^^[i(pADH-Trp-)] - (k, + kr + k,f) [3(FADH-Trp-)] + 

^fJP(FADH-*Trp-)], (5) 
d[^(FADH-Trp-)] ^ ^^j^^^.^^^.^^ _ ^^(2) ^ [i(pADH-IVp-)], (6) 
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dP(FADH-Trp')] ^ ^^p^pADH'Trp')] - (fcfj + fc^?) p(FADH-IYp-)], (7) 

=4^)[^(FADH-Trp-)] - &ADR^], (8) 

^[^^^^^ =A;sr(FADH-Trp-)] + A;Tp(FADH-Trp-)] - kred[FABR'], (9) 
"^^^^^^ ] =A;,ed[FADH'] + k^^ ([^(FADH'^Trp')] + p(FADH'*Trp')]) - 

A:(2j[FADH-]. (10) 

Here square brackets denote the concentration of transient states introduced in Fig. |2} 
[FAD], [FAD*], [i(FADH'Trp')], [^(FADH'Trp')], [i(FADH'*Trp*)], [3(FADH'*Trp')], 
[FADH**], [FADH+], [FADH*] and [FADH"] denote the concentration of FAD, FAD*, 
^[FADH* + Trp'], ^[FADH' + Trp'], ^[FADH'* + Trp'], ^[FADH'* + Trp'], FADH'*, FADH+ 
+Trp, FADH" + Trp and FADH" states, respectively. The rate constants ki, k2, k^p, kgf, 
ks, ^T, ^rel, ^rel' ^o^^ ' , Ked, ^et^ and /c^j^ in Eqs. ©-([in]) capture the different pro- 
cesses underlying cryptochrome photoactivation and relaxation shown in Fig. [21 assuming 
all processes can be described through first order kinetics. The choice of the appropriate 
rate constants is discussed in Supplementary Material and is summarized in Tab. [H 

To demonstrate that the reaction scheme in Fig. [2] is consistent with observed cryp- 



tochrome transient absorption 55] we evaluate the time development of [FAD], [FAD*], 
[i(FADH'Trp')], [^(FADH'Trp')], [i(FADH**Trp')], [3(FADH**Trp*)], [FADH'*], [FADH+], 
and [FADH"]. The evaluation can be achieved by numerical integration of the rate equa- 
tions (12|)-( !T0|) . We assume the initial condition in which holds [FAD]|j^q = cq, while all other 
intermediate states are unpopulated, cq is the initial {t = 0) concentration of cryptochrome 



in the sample before excitation by a laser pulse. According to the experimental study 55 1 
holds Co ^ 20 /iM. 



III. RESULTS 



In the following the time-dependence of the transient absorption in cryptochrome is stud- 
ied for the reaction scheme shown in Fig. [2J The calculated transient absorption spectra 
based on the scheme are then compared to available experimental data. After the postu- 
lated reaction has been validated through comparison with experimental observations, the 
influence of an external magnetic field on the signaling state of the protein is analyzed. 



11 



TABLE I. Rate constants of various kinetic processes in cryptochrome. 



Rate 


Physical process 


Char. 


Rate (s ^) 


Source 


const. 




time 






ki 


FAD excitation by blue light 


1.47 ns 


6.8 X 10^ 


estimate, see Supporting Material 


k2 


FADH* excitation by green light 


1.1 ns 


9 X 10^ 


estimate, see Supporting Material 


k 

rp 


^ [FADH' + Trp'] radical pair 










formation 


30 ps 


3.3 X 10^° 


experiment, ^3 


kef 


Singlet-H-triplet interconversion 


1 flS 


10^ 


estimate, 27 


ks 


Singlet state decay 


10 /is 


10^ 


consistent with 56 


knr 


Triplet state decay 


10 /is 


10^ 


consistent with 56 


rel 


FAD* relaxation to FAD 


80 ps 


1.25 X 10^° 


experiment, 60 


rel 


FADH** relaxation to FADH* 


80 ps 


1.25 X 10^° 


experiment, 60 




FADH"^ deprotonation 


10 ps 


lOii 


experiment, 
57, 59, 77, 79 


1.(2) 


FADH" oxidation 


4 ms 


250 


experiment, 55 


kred 


FADH* reduction 


14 ms 


70 


experiment, 55 


1.(1) 


electron transfer FADH* — > Trp* 


10 ^s 


10^ 


experiment, 56 


1.(2) 


electron transfer Trp — > FADH** 


38 ps 


2.6 X 10^° 


experiment, 60 



Characteristic time scales and rate constants for the kinetic processes in cryptochrome depicted 
in Fig. [21 



Light intensity dependence of transient spectra 

To probe the transient states in cryptochrome, a cryptochrome-containing sample from 



the retina of garden warbler had been studied (55| . The transient absorption was measured 
using the pump-probe experiment depicted in Fig. [1], where a Nd-YAG laser operating at 
355 nm excited the sample, and a 150 W tungsten lamp was used for monitoring reaction 
intermediates. To study the transient states in cryptochrome during its photoactivation cycle 
the wavelength of the probe beam was adjusted to allow absorption of either the FADH" 
and FADH* states of the flavin co-factor (absorbing at 490-550 nm), or of the FADH* state 
only (absorbing at 550-630 nm). The experimentally recorded time profiles of the transient 
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FIG. 3. Cryptochrome transient absorption spectra. Time profiles of the change in the 
transient absorption spectra from the absorbance prior laser excitation of the sample calculated 
for cryptochrome excited by laser pulses of different intensity (see Fig. [1]), for the probe-beam of 
wavelengths 550-630 nm (plot a), and 490-550 nm (plot b). The laser pulse duration chosen was 



5 ns as in experiment [55|, at 355 nm wavelength. The energies of the pulse used in the calculation 
are indicated in the inset. The wavelengths used in the calculation were chosen consistent with the 
experimental measurements of the transient absorption spectra in garden warbler cryptochrome 



55 1 . Dots correspond to the experimentally measured data recorded for the pulse energy of 5 mJ 
55|. 



absorption spectra for the two regimes are shown in Fig. [31 As discussed in the report of the 
experiment 



55| . transient absorption observed at 490-550 nm reveals a double exponential 
decay with decay times of (4± 1) ms and (14 ± 1) ms, whereas transient absorption at 550- 
630 nm shows only a single exponential decay with a decay time of (14 ± 1) ms. The decay 
times, therefore, may be attributed to FADH* (14 ms) and FADH~ (4 ms). The experimental 
error of the decay times is not provided in the original experimental papers and has been 
estimated here by digitizing the experimental data points and fitting the absorption curves 
with single and double exponential decays. 



According to the original explanation of the recorded absorption patterns, 55| the double 



exponential decay of the transient absorption spectrum was attributed to the presence of 
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the tryptophan radical in the sample for over 4 ms, and to the protein's signaling state 
associated with FADH* to be lasting over 14 ms. This explanation implies a radical pair 
lifetime in cryptochrome of over 1 ms, which is unusually long. However, a lifetime of the 
radical pair beyond 10 /is would abolish a magnetic field effect as singlet and triplet radical 
pair states ^[FADH* + Trp*] and ^[FADH* + Trp*], respectively, would then be dominated 
by electron spin relaxation that occurs faster than 10 /is j^^]. For a significant magnetic 
field effect to arise, the radical pair lifetime needs to be short compared to electron spin 
relaxation times and, in fact, needs to be comparable with hyperfine interaction-induced 
singlet-triplet interconversion times of 1 /is At present there is no physical mechanism 
known that could explain a magnetic field effect for radical pair lifetimes > 10 /is. The many 
suggestions for cryptochrome involvement in magnetoreception, for example the finding of 



98 



99| , suggests to further 



Reppert et al. for a cryptochrome knock-out mutant of fruit flies 
pursue the possibility of a cryptochrome-based biochemical compass tied likely to a radical 
pair lifetime of < 10 /is. 

who isolated 



Such short lifetime seems to contradict the finding of Brettel et al. 
cryptochrome-1 from Arabidopsis thaliana and studied the transient absorption of the pro- 
tein; analyzing the change of the transient absorption the authors suggest that a tryptophan 
radical is present in cryptochrome for over 1 ms. However, this result cannot be applied 
to the present case since Brettel et al. studied plant cryptochrome instead of avian cryp- 
tochrome. Although the photocycle of the members of the cryptochrome family should 
exhibit similarity, experiments sh ow g reat variation between cryptochromes extracted from 



different organisms 77|, l80|, l83|, |84, l86|, |87| . Furthermore, the change in the transient absorp- 



tion spectra attributed by Brettel et al. to the presence of the tryptophan radical may stem 
from another cryptochrome component with similar absorption properties. To explain an 
error due to such component we note that cryptochrome transient absorption relative to the 
absorbance prior to laser excitation can be calculated using Eq. ([1]), where the summation 
is performed over the FAD, FADH~ and FADH* states for a 490-550 nm probe-beam and 
includes only one term corresponding to the FADH' state in case of a 550-630 nm probe- 
beam. In the suggested model the transient ^[FADH* + Trp*], ^[FADH* + Trp*] states have 
no noticeable influence on the absorption spectrum since the lifetime of these radical pair 
states is assumed to be several orders of magnitude shorter than the lifetime of the FADH" 
and FADH* states. If the lifetime of the radical pair state is on the order of milliseconds, as 
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suggested in (1551 ). then the absorbance of the tryptophan radical becomes important. 

Figure [3] shows the calculated time- dependence of the transient absorption in cryp- 
tochrome, obtained from solving rate equations fj2|)-f fT0|) . for wavelengths at 490-550 nm 
and 550-630 nm. The figure shows that the energy of the laser pulse impacts the absorp- 
tion properties of the protein. The dependence on the laser-pulse energy is due to the rate 
constant /c2 which is proportional to the laser pulse power (see Supplementary Material for 
more details). The transient absorption spectrum calculated for a pulse power of 1 MW, 
i.e., corresponding to the experimental 5 ns pulse carrying 5 mJ energy, can be compared to 
the experimental result reported earlier 55| . As can be seen in Fig. [3l the calculated tran- 
sient absorption perfectly matches the recorded data for a 490-550 nm and a 550-630 nm 
probe-beam. 

Increasing the energy of the laser pump-pulse enhances the rate constant ^2 leading 
to faster population of the FADH~ state (see Fig. |2]). This enhancement can be verified 
experimentally since the FADH" and FADH' states in cryptochrome absorb light differently. 
Indeed, exciting the cryptochrome-containing sample with a pulse of higher intensity should 
lead to diminished absorption at 550-630 nm and a single exponential decay at 490-550 nm 
with decay time of 4 ms. Figure |3] illustrates how the time-dependence of the transient 
absorption behaves once the intensity of the pump pulse is increased. 

A decrease of the pump-pulse power leads to suppression of the FADH* — )■ FADH** 
FADH~ channel and extends the lifetime of the [FADH* + Trp*] radical pair. Figure |3] shows 
the dependence of the transient absorption calculated for the laser pulse energy of 1 mJ act- 
ing over 5 ns. For such weak and short pulse the transition FADH* — )■ FADH** — )■ FADH" 
is significantly suppressed as the redistribution of the FADH* concentration towards the 
FADH" state is quenched. Figure [3] shows that for lower light intensities the absolute mag- 
nitude of the absorption at 550-630 nm increases, because the concentration of the FADH* 
states increases. Figure |3] presents experimental data only for laser pulse energy of 5 mJ; 
to establish unambiguity of the suggested model it is necessary to perform measurement at 
laser pulse energies of 1, 4, 7 and 10 mJ. 

It is worth noting that the photoexcitation rate constants ki and k2 depend not only 
on the laser-pulse energy, but also critically on the wavelength of the pulse and on pulse 
duration. Figure H] presents the time evolution of the concentration of the transient states 
in cryptochrome calculated from Eqs. fl2|)-( fT0|) for the photoexcitation regimes considered 
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FIG. 4. Population of transient states in cryptochrome. Time evolution of concentration of 
the transient states in cryptochrome calculated from equations (|2])- (ll0p for cryptochrome photoex- 
citation b y a laser pulse of 5 ns duration, 5 mJ energy and 355 nm wavelength (plot a, consistent 



with Ref. 



55l ). and a laser pulse of 6 ns duration, 4 mJ energy and 460 nm wavelength (plot b. 



consistent with Ref. 



56l ) . Grayed areas indicate time duration of laser pulse. 



in one study 55| (plot a) and in another 56|] (plot b). The difference in the pump- pulse 
leads to a significant change in the population of transient states. Thus, for a laser pulse 
of 5 ns duration, of 5 mJ energy, and 355 nm wavelength (see Fig. HH) the population of 
the FADH^ state after 5 ns is about 5 times higher than the population of the radical pair 
^[FADH* + Trp*] state. The early population of the FADH~ state readily explains both the 
double-exponential decay of the transient absorption spectrum at 490-550 nm and a single 



exponential decay at 550-630 nm as reported earlier 



55|. 



The absorption of the semiquinone state of the flavin radical (FADH*) at 460 nm is 
significantly weaker in comparison to absorption at 355 nm (see Fig. [7] in Supplementary 
Material), and, therefore, the photoexcitation rate constant /c2 at 460 nm is smaller than at 
355 nm. Figure Hb shows that for a laser pulse of 6 ns duration, 4 mJ energy, and 460 nm 
wavelength the population of the FADH" state after pulse excitation is smaller than the 
population of the radical pair ^[FADH* + Trp*] state, arguing that the radical pair state can 
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be better resolved at 460 nm than at 355 nm wavelength. 



Influence of magnetic field as seen by transient spectra 



Comparison of measured and calculated kinetics in cryptochrome reveals good agreement 
and, hence, supports the reaction scheme in Fig. [2l Therefore, the postulated reaction 
scheme will be employed further to consider magnetic field effects on cryptochrome activation 
and deactivation. 

The influence of an external magnetic field on cryptochrome photoactivation can be 
probed through the effect of an applied magnetic field on the transient absorption spectrum. 
The magnetic field influences the singlet-triplet interconversion process in cryptochrome (see 
Fig. [2]), described here extremely schematically through the first-order rate constant kgf- An 



applied magnetic field leads to a relative change of the rate constant kgf 24 
brings about an altered transient absorption. The change can be written 



25 



46[ which 



AA{a,t) = A{ksf,t)-A{ak,f,t), (11) 

where kg/ = 10^ s~^ is the reference value of the singlet-triplet interconversion rate constant 
(see Tab. [11) and a denotes the relative change in kgf due to an altered magnetic field; A is 
defined in Eq. ([1]). 

The a-value depends on several factors, such as the hyperfine interaction in the radical 
pair, exchange and dipole-dipole interactions, and Zeeman interaction 24j. |25| . An acc urate 



calculation of a requires the solution of a stochastic Liouville equation 
consider a- values ranging from a = 1.2 [k^f = 1.2 x 10^ s~^) to a = 



24 



25 



46 



100|. We 

10^ s-M, 



10 {ksf 

thereby accounting for realistic singlet-triplet interconversion rates expected in cryptochrome 



when the field reorients 



24 



25 



27| 



Deviation of a from a = 1 can arise either through a change of magnetic field orientation 
or through changing field strength. The latter change can arise due to geographic variation 
of the geomagnetic field strength; this change is small percentage- wise, but is employed 
by birds apparently for a so-called magnetic map sense 10ll - ll04| . In this case detection 



is like 
beak 



y achieve d by an alternative receptor, an iron-mineral-based magnetoreceptor in the 



104 



106|. The geographic variation of the magnitude of the geomagnetic field is 



expected to have a negligible impact on the photo-reaction kinetics of cryptochrome and. 



17 




FIG. 5. Magnetic field effect in cryptochrome. Change of the transient absorption spectra 
in cryptochrome due to a change in kgf caused by an external magnetic field, calculated for the 
550-630 nm probe-beam, a is defined in the text. Initial photoexcitation is due to a laser pulse 
of 5 ns duration, 355 nm wavelength and 5 mJ energy (plot a) and 2.5 mJ energy (plot b). The 
curves in (a) and (b) reflect change in transient absorption due to an increased value of the kgf 
relative to the reference value 10® s~^ (see Tab. [I]). 

therefore, is not considered here; only a change, i.e., a ^ 1, due to reorientation in the 
Earth' field is considered. 

For cryptochrome to exhibit a response to the change of the magnetic field direction 
the protein needs be constrained orientation-wise within the organism. We have recently 
demonstrated that only one of three rotational degrees of freedom of cryptochrome need 



to be constrained to endow a bird with the magnetic compass [26|]. Such constraint can 
be realized if cryptochromes are embedded in a cell membrane. The outer segment of the 
photoreceptor cells is an ideal structure to constrain the proteins as originally suggested in 



theory 26|. An experimental verification of cryptochromes localization in the outer segment 



67|. 



of the UV photoreceptor cells was published recently 

Figure [5] shows the time- dependence of AA for various a- values. The calculations were 
performed for cryptochrome photoexcitation by a laser pulse of 5 ns duration operating at 
355 nm and having the energy of 5 mJ (see Fig. [5^) and 2.5 mJ for which the life-time of 
the [FADH* -|- Trp*] radical pair is expected to be longer (see Fig. [5]d). The relative change 
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of ksf due to an applied magnetic field, ^ ma y vary within an order of magnitude as stated 
above and in Supplementary Material 



24 



27|. 



Figure [5] shows that decreasing the intensity of the pump-pulse leads to an enhancement 
of the magnetic field effect. The enhancement is due to the excitation channel FADH* — j- 
FADH** — 7- FADH~ in the cryptochrome photoactivation reaction (see Fig. [2]) which becomes 
suppressed for the 2.5 mJ energy laser-pulse, such that the lifetime of the [FADH* + Trp*] 
radical pair increases. Increasing ksf {a > 1) leads to faster singlet-triplet interconversion 
in the radical pair. If this interconversion becomes significantly faster than singlet and/or 
triplet decay, governed by rate constants ks and /ct (see Fig. [2]), the magnetic field effect in 
cryptochrome becomes negligible, as clearly seen in Fig. [5]d, where the relative difference in 
AA calculated for a = 5 and a = 10 is small. A similar behavior is expected for higher pump- 
pulse intensities (see Fig. [5^) at higher values of a, because in this case the [FADH* + Trp*] 
radical pair decays predominantly through the FADH* — )■ FADH** — j- FADH" channel, 
governed by the rate constant k2, which is significantly larger than ks and /ct- 

The dependencies shown in Fig. |5] can also be measured experimentally and are important 
for resolving the physical mechanism of the magnetic field effect in cryptochrome. Similar 



59| and should be repeated for 



measurements were performed earlier for DNA photolyase 
cryptochrome. The time- dependence of the transient absorption change due to a magnetic 
field allows one not only to judge whether the protein exhibits magnetoreceptive proper- 
ties, but can also be used to reveal the impact of the magnetic field on the singlet-triplet 
interconversion rate constant kgf. 

Predicted absolute magnetic field effects in cryptochrome are rather small, i.e., the ab- 
sorption changes only by about (0.1-0.3) x 10~^ for a 5 mJ pulse (see Fig. [5|i) and by about 
(0.2-0.6) X 10"^ for a 2.5 mJ pulse (see Fig. |5b). Absolute changes of the transient absorp- 
tion due to a 39 mT magnetic field in DNA photolyase were recorded to be of the order 
of ~0.001 59|, e.g., larger, but still very small. The relative magnetic field effect for DNA 



photolyase can thus be estimated to be ~1 %. The difference in the magnetic field effect 
recorded for DNA photolyase and predicted now for cryptochrome is due to the concentra- 
tion of DNA photolyase used in experiment 59|], which was 0.2 mM, whereas in the present 
study the assumed concentration of cryptochrome is only 20 yuM, following an earlier study 



55(1 . Increasing the concentration leads to an enhancement of the absolute magnetic field 
effect. 
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FIG. 6. Magnetic field effect and probe concentration. The maximal change of the transient 
absorption spectra in cryptochrome calculated due to change in kgf at different concentrations of 
protein, calculated for a 550-630 nm probe-beam. Initial photoexcitation is due to a laser pulse of 
5 ns duration, 355 nm wavelength, 5 mJ energy (plot a) or 2.5 mJ energy (plot b). The curves 
in (a) and (b) show transient absorption due to an increased k^f relative to the reference value of 
10^ s~^ (see Tab. H]), a is defined in the text. 

In an animal the ambient light continuously irradiates cryptochromes in the retina and 
photo-excitation of FADH* induced through strong irradiation does not arise because the 
rate constant k2 becomes small as compared to the value in the transient absorption ex- 
periment. In vivo conditions lead then to an increase of the [FADH* + Trp*] radical pair 
lifetime, increasing also the effect of an external magnetic field on cryptochrome signaling 
associated with FADH*. 

Figure |6] shows the maximal change of the transient absorption due to a magnetic field 
effect calculated at different concentrations of the sample, and at different intensities of the 
pump-pulse (Fig. [6^ for 5 mJ over 5 ns and Fig. [6]d for 2.5 mJ over 5 ns). From the analysis 
in Fig. |6] follows that the magnitude of the maximal change of the transient absorption in 
cryptochrome increases linearly with the concentration of the sample. Thus, increasing the 
concentration to 200 /iM could lead to a magnetic field effect of about 0.0007 (see Fig. Eb), 
i.e. close to the value detected for DNA photolyase jsol. It is also worth noting that a more 
realistic account of the singlet-triplet interconversion predicts a relative magnetic field effect 
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in cryptochrome of about 5-10 % in a magnetic field comparable to the geomagnetic field 



24 



25| . which suggests that a is large, i.e. a ^ 1, and that the reference value of kgj is 



< 10^ s-i. 

The trivial concentration dependence of the maximal magnetic field effect in cryp- 
tochrome shown in Fig. [6] carries important information about the rate constants in cryp- 
tochrome photoactivation reaction. In particular, the slopes of the linear dependencies 
shown in Fig. [6] depend on these rate constants which govern the lifetime of transient states 
in the photoactivation reaction. The dependencies shown in Fig. [6] are convenient for the 
experiment, because the concentration of the cryptochrome-containing sample can be easily 
varied, and the transient absorption of the sample can then be recorded. Thus, by com- 
paring results of the measurement with predictions of the kinetic model suggested in this 
paper it is feasible to justify rate constants in the photoactivation reaction, and predict the 
lifetime of certain intermediate states in cryptochrome. 

The amount of cryptochrome in the bird's retina, presently unknown, is critical for the 
magnitude of the magnetic field effect. For a maximal signal a cryptochrome concentration 
larger than 20 /iM is necessary. Figures [5] and [6] allow one to draw another important 
conclusion from our analysis, namely, that decreasing the intensity of the pump-pulse leads 
to an enhanced of the magnetic field effect. Interestingly, many migratory birds favor to fly 
at night time, i.e., when the intensity of the ambient light is low 



IV. CONCLUSION 



Earlier investigations 



23-26 



45l | demonstrated how a biochemical magnetic compass in 



birds can be realized through light-induced electron transfer reactions in cryptochrome. 
However, the biophysical mechanism of cryptochrome activation and signaling still remains 
elusive. Our study provides an important step towards resolving the mechanism of cryp- 
tochrome magnetoreception. The suggested analysis can provide a proof-of-principle whether 
cryptochrome is indeed the magnetoreceptor proteins in birds, comparing computed kinetic 
data with data from transient absorption spectroscopy, as successfully applied to DNA pho- 



tolyase 



59|. 



The results of the present calculations reveal that magnetic field effects become only dis- 
cernable if the concentration of the protein in the probe exceeds a critical value. Therefore, 
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cryptochrome should be measured at different concentrations to ascertain that the minimal 
amount of protein is available. Transient absorption of cryptochromes at low concentra- 
;ion can be studied with the use of the so-called cavity-enhanced spectroscopy technique 
. This technique traps the probe light beam in an optical cavity, containing the sample, 



thereby transmitting the beam through the sample multiple times and increasing the trans- 
mittance length which in turn enhances the total absorption in the sample. The method 



permits protein concentrations below ~1 fiM [83[, resolving nevertheless small changes in 
the transient absorption, as they arise due to an external magnetic field. A systematic 
study of the dependence of the absolute absorption of a cryptochrome-containing sample on 
cryptochrome concentration is important for revealing the lifetime of reaction intermediates 
through direct comparison of experimentally recorded absorption with the prediction of the 
kinetic formalism suggested in this paper. 

The methodology suggested here illustrates how laser power, laser wavelength, laser pulse 
duration, probe light wavelength, and protein concentration, impact measured results. In 
particular, we identify which measurements pinpoint the magnetoreceptive properties of 
cryptochrome. 

The performed study assumes for its analysis reaction rate constants, which were ex- 
tracted from experiment or could be estimated otherwise. The uncertainty of most of the 
rate constants, is rather large. By measuring the transient absorption spectra in cryp- 
tochrome at different conditions, i.e., for different probe-beam intensities, different applied 
magnetic fields, different wavelengths of the pulse and probe beams, etc, and by comparing 
recorded profiles with calculated ones one can narrow the range of consistent rate constants 
in the proposed reaction scheme. 

Magnetoreception in birds has been found to work for light with wavelengths of up to 565- 



567 nm 
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recep 

iqtI. 



107| . Thus, it is important to continue studies of cryptochrome photoactivation 



in v itro by maintaining the conditions close to those used in the in vivo experiments 
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107j . An important question is how the magnetic directional information is disrupted under 



monochromatic light of an increased intensity. It has been d emonstrat ed that different 



light intensities can disorient birds in the laboratory 18 



20 



108 



109| . and, therefore. 



light intensity may affect the photocycle in cryptochrome. Studying transient absorption 
of cryptochrome in electromagnetic fields of certain frequencies, e.g., Larmor frequency of 



the electron 



42 



43 



96j . is another step to link in vitro studies to behavioral investigations 
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42|,|43|. 

Science is still a long way from an explanation of avian magnetoreception. The best that 
may be said of our understanding today is that birds indeed perceive and use magnetic field 
information, and that their responses to magnetic fields under different conditions - light 
intensity and color, magnetic field strength and presence of oscillating fields - are suggesting 
a complex sensory system of multiple receptors. Even if cryptochrome is confirmed as a 
magnetoreceptor, it remains for biologists to determine how its magnetic field modulated 
signaling enters into a bird's sensory perception and ultimately into its orientation behavior. 



V. SUPPORTING MATERIAL 

In Supporting Material we present details on flavin absorption spectra and discuss the 
choice of the rate constants in the studied reaction in cryptochrome. 



A. Absorption spectrum 



The efficiency of light absorption at wa velen gth A by an absorbing medium is characterized 
through the absorbance A{X) defined as jllO | 

/o(A) 



AiX) = log 



/(A) 



(12) 



where Io{X) and /(A) are the light intensities of the beam entering and leaving the absorbing 
medium, respectively. According to the reaction scheme in Fig. |2l cryptochrome can occupy 
several states, which are expected to absorb light differently. Therefore, according to the 



110|], the absorbance of the cryptochrome-containing sample is 



Beer-Lambert-Bouguer law 
given by Eq. (P). 

The molar extinction coefficient, Ei, in Eq. ([T]) of a light-absorbing component in the 
system is directly related to the absorption cross section, (Xj, which characterizes the photon- 
capture area of a molecule [110 1 



ai{X) = lOOOln(lO)- 



Na 



3.82 X 10-^^ei{X). 



(13) 



Here Na is Avogadro's number; the absorption cross section in Eq. (fT3|) is measured in units 
of cm^. 
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B. Rate constants 



The rate constants in Eqs. (!2|)- (lT0l) determine the time evolution of intermediate states 
in cryptochrome. Many of the rate constants are available from experiments performed on 
cryptochrome and cryptochrome-like proteins from various species, such as garden warbler 



55 1 . Drosophila melanogaster 79|, |80[, Arabidopsis thaliana 77|, l78l|, and Homo sapiens 



56|,Q- 

Many studies were also done for DNA photolyase, a protein structurally similar 



the rate constants can be independently estimated 



2i 



25|. The adopted rate constants are summarized 



to cryptochrome 57H6ll. 1631 ] . Some o 
from fundamental physical principles 
in Tab. [B 

Flavin photoexcitation. The rate constants ki and k2 represent the rate of photoexcita- 
tion of the flavin cofactor from its fully oxidized FAD state and from the semiquinone FADH* 
state to the excited FAD* and FADH** states, respectively, (see Fig. [2]). The FAD -> FAD* 
and FADH* — )■ FADH** transitions are induced by a laser pulse (see Fig. [1]) and arise only 
during the pulse duration r. The rate constants ki and k2 depend on the laser power and 
can be estimated as 



fc,. = a/^(l-e(t-r)). 



(14) 



where P = E/t is the power of the laser pulse (with E being the energy of the pulse and r 
the pulse duration), 5*0 is the cross section area of the light beam hitting the sample, x < 1 
defines the fraction of power deposited at the sample, a is the absorption cross section 
defined in Eq. fll3p . Eph = hc/X is the energy of a single photon (with h being the Planck 
constant and c the speed of light), 6(x) is the Heaviside step-function which limits the 
photoexcitation of FAD and FADH* to the period of the laser pulse duration. Here we do 
not consider the periodicity of the laser pulses as the time interval between two successive 
pulses is significantly longer than the typical reaction times involved in the scheme shown 
in Fig. [2] (the pulse frequency used in the observation [55| was 10 Hz). 

Assuming a Gaussian radial profile of the laser beam and substituting Eq. ( fT3l) into 
Eq. ([T^ one obtains 



612.686 X 



PXeiX) 

i?2 



exp 
exp 



(l-e(t-r)) 



(15) 
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where R is the radius of the beam at the sample measured in /im, Rq is the radius of the 
output beam from the laser measured in fim and u is the radius at which the laser field 
amplitude drops to 1/e. e in Eq. ( IT5l) is measured in L ■ mol~^cm~^, A is measured in nm 
and P is measured in Watt. 

In the measurements jsSi the Rq value was Rq = 3000 fim and the size of the sample 
used was likely smaller allowing one to assume R = 15 00 /xm. W ith u = 1000 fim, a 



typical value for the amplitude fall-off of the laser beam [ill 



112) ■ Eq. ([T5]) can be used 



to estimate the photoexcitation rate constants ki and The rate constants fci ^-nd are 



determined by the laser wavelength, 355 nm, and the beam power, 10^ W 



55| . Figure [7] 



shows that the extinction coefficient of the oxidized flavin (FAD) at A = 355 nm is about 
7900 L ■ mol-^cm-i 



55 



58 



resulting in ki = 6.8 x 10^ s~^. The absorption (extinction) 
spectra in Fig. [7] were recorded for the three redox states of FAD and normalized to the 
absorption of FAD at 450 nm (£450 = 11.3 x 10"^ L ■ mol~^cm~^) (solid line) 55[. We note that 



the wavelength dependence of the extinction coefficient in garden warbler 



and human 



56 



55j . Drosophila, 



80j cryptochromes maintains the general features of the absorption profile 
shown in Fig. [71 Similarly, Fig. [7] shows that the extinction coefficient of semiquinone 



(FADH*) flavin at A = 355 nm is about 10400 L ■ mol ^cm ^ (solid line) 55|, |58[, resulting 



in /c2 = 9 X 10^ s ^. The rate constants ki and /c2 apply only during the laser pulse duration 
time of r = 5 ns. 



Flavin excited states relaxation. The rate constants and fc^g^ describe the FAD* 
FAD and FADH** — )■ FADH* relaxation processes, respectively (see Fig. [2]). These relaxation 
processes have not been very well documented in cryptochrome, but the FADH** FADH* 



,(2) 



transition was studied in DNA photolyase 



60|^ protein homologous to cryptochrome 



113] . According to experimental measurement 60| the lifetime for the relaxation of FADH** 



to the ground state is 80 ps, leading to the value fc^l = 1-25 x 10^° s~^ The FAD* FAD 



transition is expected to occur on a similar timescale and, accordingly, we assume k. 



(1) 

rel 



1.25 X 10^" s~\ 

Radical pair formation. The rate constant krp describes the ^ [FADH* + Trp*] radical 
pair formation process (see Fig. |2]). The characteristic time for this process is 30 ps, as 
confirmed by using ultrafast pump-probe spectroscopy in the near-infrared spectral region 
in DNA photolyase 63|, leading to krp = 3.3 x 10^'^ s~^. Although krp has not been clearly 
resolved for cryptochrome, we assume the rate constant from DNA photolyase to be of the 



25 



same order of magnitude in all photolyase/cryptochrome-like proteins 
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56 



59|. 



Flavin deprotonation. The rate constant ki^J describes the deprotonation process of 



FADH+, as denoted in Fig El Since the FADH 
observed in cryptochrome, and/or photolyase 



state o: 



55-5 



3 



the flavin cofactor has never been 



59-61 



63, 



77-79 



114l |. the character- 



istic time of FADH"*" deprotonation is expected to be on the order of few picoseconds, which 
is beyond experimental resolution. Thus, we assume the value k^J = 1/10 ps = 10^^ s~^. 

Electron transfer involving flavin radical. The rate constants /CgJ'* describes electron 
transfer from the FADH* radical to the Trp* radical (see Fig. |2]). k^^^ is expected to have an 
approximate value of 1/10 fis = 10^ s~^, since the radical pair in cryptochrome is assumed 



to have a lifetime of 6 fis |56|; measurements were performed using transient EPR spec- 
troscopy, with the system optically excited by a Nd:YAG laser (Spectra Physics GCR-11) 
pumping an optical parametric oscillator (Opta BBO-355-vis/IR, Opta GmbH, Bensheim, 
Germany) tuned to a wavelength of 460 nm (pulse width 6 ns; pulse energy 4 mJ). At this 
particular wavelength the extinction coefficient of the semiquinone (FADH*) flavin is about 
3600 L ■ mol~^cm~^ (see blue solid line in Fig. [7]). The chosen laser pulse power leads to the 
photoexcitation rate constant k2 = 2.6 x 10® s~^ under the assumption that the geometrical 
characteristics of the laser beams in two experiments 55|, l56|| are identical. Since the fast 
photoexcitation decay channel of the radical pair is only possible during the period of the 
laser pulse, the fourfold decrease of the photoexcitation rate constant would lead to a sig- 
nificant increase of the lifetime of the radical pair state allowing its detection in the EPR 
measurements. 

(2) 

The rate constant k^^ describes electron transfer from the Trp* radical to the excited 
FADH** radical. The electron transfer rate constant for this process was measured in DNA 



ohotolyase using time-resolved absorption spectroscopy and found to be k^'' = 2.6 x 10^*^ s~^ 
60(1 . Similar values for the electron transfer rate constants can also be estimated from Marcus 
theory of electron transfer 24 1. 

S ^ T interconversion. The rate constant ksf describes the singlet-triplet interconver- 
sion process in cryptochrome (see Fig. [2]) as a first order reaction process. We employ such 
process as a very rough model for the actual quantum mechanical spin precession process, 
since only time scale and yield of singlet-triplet interconversion matter. The ksf rate constant 
depends on several factors, such as the hyperfine interaction in the radical pair, exchange 
and dipole-dipole interaction between the radical pair partners, and the external magnetic 



26 



field. To calculate the singlet-triplet kinetics one nee ds to solve the stochastic Liouville 



equation for the radical pair in the system 



25 



46 



100|. In the present study we cast 
It has been demon- 



the transition process into a single rate constant as done earlier 

n n □ 

strated [24, |25|, |27| that for a generic radical pair holds ksf = 10^ — 10^ s ^. kgf is magnetic 
field dependent and, thereby, responsible for the magnetic field effect in cryptochrome. We 
demonstrate the feasibility of a magnetic field effect in cryptochrome by varying the value 
of ksf. If cryptochrome is the primary magnetoreceptor protein in birds and other animals, 
it is natural to assume that nature has designed it in such a way that the external geomag- 
netic field produces a significant effect. The time needed for a significant transformation 
of a singlet state ^[FADH* + Trp*] into a triplet state ^[FADH* + Trp*] and vice versa in a 
0.5 G magnetic field is typically ~ 700 ns 27|; therefore we assume kg/ = 10^ s~^ for the 
singlet-triplet interconversion rate constant. 



Singlet and Triplet decay kinetics. The rate constants kg and k^ describe the singlet 
and the triplet decays of the [FADH* + Trp*] radical pair (see Fig. |2]). Together with the 
fcgj'' electron transfer rate constant they define the lifetime of the radical pair. According 



to the experiment, the lifetime of the radical pair in cryptochrome is > 6 /is 



; therefore, 



we assume ks = kx = 1/10 /is = 10^ s~^, i.e. we take the lower bound as our estimate of 
the radical pair lifetime. For the sake of simplicity we also assume spin-independent decay 
kinetics of the radical pair, i.e. ks = /cr- 



Dark reaction kinetics. The rate constants kred and kox are associated with the two- 
step cryptochrome dark reaction (see Fig. [2]). Both stages involved are expected to be 
fairly long-lived, with a lifetime on the order of milliseconds, and should be resolved in 



transient absorption spectra reported in the experiment [55|. According to the experiment, 
two transient states in cryptochrome with lifetimes of 4 ms and 14 ms were detected, kred 
describes the reduction process of the semiquinone FADH* radical possibly involving the O'" 
radical, as suggested earlier 25| . It is natural to expect that the lifetime of the signalling 
state in cryptochrome is maximal, leading to the values k^ed = 1/14 ms ~ 70 s~^ and 
kf] = 1/4 ms = 250 s~^. 
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C. Flavin absorption spectra 



Optical absorbance of cryptochrome is dominated by the absorbance of the isoalloxazine 



55 



moiety in FAD 
states of FAD: the fn 
oxidized FAD 



56| . Figure [7] shows absorption spectra recorded for the three redox 



56 



89 



ly reduced FADH , one-electron oxidized FADH* and the two-electron 



97|. 




FIG. 7. Flavin absorption spectra. Wavelength dependence of the molar absorption of the 



^ADH" 



55| and 



flavin chromophore in its three redox states FAD (red line), FADH* (blue line), and 
(green line). Vertical dashed lines set the limits for the wavelengths used in experiment 
in the calculations reported here. The spectra of the flavin chromophore shown as a solid line 
have been adopted from a textbook [971], while the spectra of different redox forms of selected 



cryptochrome/photolyase shown as a dashed line are digitized from the measured pattern 



The calculated absorption spectrum in Fig. |3] depends on the value of the extinc- 
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tion coefficients, Si, of the cryptochrome transient states (see Eq. ([T])), which are wave- 
length dependent. The extinction coefficients for FAD, FADH' and FADH~ were chosen 
as 300 L ■ mol~^cm~^, 2695 L ■ mol~^cm~^ and 690 L ■ mol~'^cm~^, respectively, for 490- 
550 nm, and the extinction coefficient for FADH* at 550-630 nm was assumed to be 
3843 L ■ mol~^cm~^. Figure [7] illustrates that the chosen values are consistent with the 
experimentally recorded absorption profiles for the different redox states of the flavin moi- 
ety. 

The only noticeable difference between calculated and observed spectra is the increased 
value of the extinction coefficient for the fully reduced FADH^ state of the flavin cofactor, 
which in the calculation is slightly larger than the extinction coefficient for the fully oxidized 
FAD state. However, both extinction coefficients are expected to be small (see Fig. [7]), as 
they are taken from the far edge of the absorption spectrum. Therefore, the inaccuracy in 
their value is high, but inconsequential. To illustrate the uncertainty of flavin extinction 
at 490-550 nm in Fig. [7] we show two sets of absorption profiles for the different redox 
states of the flavin cofactor recorded for an isolated flavin chromophore (solid line) 
and for the selected cryptochrome/photolyase proteins 89| (dashed line). Although in both 
systems the absorption spectra contain similar features, one notes that the spectra for the 
oxidized (FAD) and semiquinone (FADH*) states in the isolated flavin chromophore are red- 
shifted with respect to the spectra recorded for cryptochrome/photolyase proteins, thereby 
exhibiting significantly different absorption at 490-550 nm. Another reason for the increased 
absorption of the FADH* state in our model may be that in the experiment 5^ FADH" is 
actually substituted by the anionic FAD*~ radical state which exhibits a somewhat higher 



absorption at 490-550 nm 



891]. This explanation is worth studying in a systematic fashion. 



Transient absorption is a powerful technique for studying intermediate states in chemical 
reactions, but the method often delivers data which cannot be unequivocally interpreted 
since the absorption patterns for different molecules and their various states often overlap. 
The impact of simultaneous absorption of several transient states at the same wavelength 
on the total absorption of the sample had to be addressed already three decades ago for 
pyrene/DMDMA complex [s^], where a theoretical approach, similar to the method which we 
now use for cryptochrome transient absorption calculation, was successfully employed in the 
ffist quantitative demonstration of a radical pair magnetic field effect. Large biomolecules 
usually have many constituents which respond to the typical wavelengths used by probe 
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light in the transient absorption experiment. Recent studies on transient absorption of 
cryptochrome, and cryptochrome-hke proteins, often dehver a single explanation for the 
measured data and do not discuss how the result would change if other possible transient 
states of the studied proteins were taken into account. The only way to unequivocally 
interpret the results of measurements in the present case is through measurements at different 
conditions that likely impact the reaction kinetics and shift the equilibrium toward different 
transient states. 

ACKNOWLEDGMENTS 

The authors thank Paul Galland (University of Marburg) for providing the necessary data 
on flavin absorption, Robert Bittl (Freie Universitat Berlin) for stimulating discussions which 
helped us understand the experimental aspects of the transient absorption spectroscopy, 
Peter Hore and colleagues (University of Oxford) for providing the necessary information 
on the experimental setup used in cryptochrome transient absorption experiments and Vita 
Solovyeva for inspiring remarks. This work has been supported by National Science Founda- 
tion grants NSF MCB-0744057 and NSF PHY0822613 and by National Institutes of Health 
grant P41-RR005969. K.S. thanks the Alexander von Humboldt Foundation for support 
and I.S. acknowledges support Beckman Fellow. 



[1] H. Mouritsen and T. Ritz, Curr. Opin. Neurobiol. 15, 406 (2005). 

[2] R. Wiltschko and W. Wiltschko, BioEssays 28, 157 (2006). 

[3] I. A. Solov'yov, K. Schulten, and W. Greiner, Physik Journal 9, 23 (2010). 

[4] W. Wiltschko and R. Wiltschko, Science 176, 62 (1972). 

[5] F. Merkel and W. Wiltschko, Vogelwarte 23, 71 (1965). 

[6] R. Freire, U. Munro, L. J. Rogers, S. Sagasser, R. Wiltschko, and W. Wiltschko, Animal 

Cognition 11, 547 (2008). 
[7] J. Voss, N. Keary, and H. Bischof, Neuroreport 18, 1053 (2007). 
[8] W. Keeton, Proc. Natl. Acad. Sci. USA 68, 102 (1971). 
[9] W. Wiltschko and R. Wiltschko, Nature 291, 433 (1981). 



30 



[10] W. Wiltschko and R. Wiltschko, Naturwissenschaften 85, 164 (1998). 

[11] G. Fleissner, B. Stahl, P. Thalau, G. Falkenberg, and G. Fleissner, Naturwissenschaften 94, 
631 (2007). 

[12] I. A. Solov'yov and W. Greiner, Biophys. J. 93, 1493 (2007). 
[13] I. A. Solov'yov and W. Greiner, Eur. Phys. J. D 51, 161 (2009). 
[14] I. A. Solov'yov and W. Greiner, Phys. Rev. E 80, 041919 (2009). 

[15] J. L. Kirschvink, M. M. Walker, and C. E. Diebel, Curr. Opin. Neurobiol. 11, 462 (2001). 
[16] G. Falkenberg, G. Fleissner, K. Schuchardt, M. Kuehbacher, P. Thalau, H. Mouritsen, 

D. Heyers, G. Wellenreuther, and G. Fleissner, PLoS ONE 5, e9231 (2010). 
[17] W. Wiltschko, U. Munro, H. Ford, and R. Wiltschko, Nature 364, 525 (1993). 
[18] W. Wiltschko and R. Wiltschko, J. Exp. Biol. 204, 3295 (2001). 

[19] W. Wiltschko, A. Moller, M. Gesson, C. Noll, and R. Wiltschko, J. Exp. Biol. 207, 1193 
(2004). 

[20] R. Muheim, J. Backman, and S. Akesson, J. Exp. Biol. 205, 3845 (2002). 
[21] R. Wiltschko, U. Munro, H. Ford, K. Stapput, and W. Wiltschko, J. Exp. Biol. 211, 3344 
(2008). 

[22] K. Schulten, C. E. Swenberg, and A. Weller, Z. Phys. Chem. NFlll, 1 (1978). 

[23] T. Ritz, S. Adem, and K. Schulten, Biophys. J. 78, 707 (2000). 

[24] I. A. Solov'yov, D. Chandler, and K. Schulten, Biophys. J. 92, 2711 (2007). 

[25] I. A. Solov'yov and K. Schulten, Biophys. J. 96, 4804 (2009). 

[26] I. A. Solov'yov, H. Mouritsen, and K. Schulten, Biophys. J. 99, 40 (2010). 

[27] C. T. Rodgers and P. J. Hore, Proc. Natl. Acad. Sci. USA 106, 353 (2009). 

[28] K. Salikhov, Y. Molin, R. Sagdeev, and A. Buchanenko, Spin polarization and magnetic field 

effects in radical reactions, vol. 22 (Budapest: Elsevier, 1984). 
[29] U. Steiner and T. Ulrich, Chem. Rev. 89, 51 (1989). 
[30] R. Beason and P. Semm, J. Exp. Biol. 199, 1241 (1996). 

[31] W. Wiltschko, U. Munro, H. Ford, and R. Wiltschko, Proc. R. Soc. Lond. B. (Biol. Sci.) 
276, 2227 (2009). 

[32] W. Wiltschko, U. Munro, R. Beason, H. Ford, and R. Wiltschko, Experientia 50, 679 (1994). 
[33] M. Zapka, D. Heyers, C. M. Hein, S. Engels, N.-L. Schneider, J. Hans, S. Weiler, D. Dreyer, 
D. Kishkinev, J. M. Wild, et al.. Nature 461, 1274 (2009). 



31 



[34] K. Schulten and A. Weller, Biophys. J. 24, 295 (1978). 

[35] K. Schulten and P. G. Wolynes, J. Chem. Phys. 68, 3292 (1978). 

[36] K. Schulten, in Biological Effects of Static and Extremely Low Frequency Magnetic Fields, 
edited by J. H. Bernhard (MMV Medizin Verlag, Munich, 1986), pp. 133-140. 

[37] K. Schulten, H. Staerk, A. Weller, H.-J. Werner, and B. Nickel, Z. Phys. Chem. NFlOl, 371 
(1976). 

[38] K. Schulten, J. Chem. Phys. 80, 3668 (1984). 

[39] H.-J. Werner, K. Schulten, and A. Weller, Biochim. Biophys. Acta 502, 255 (1978). 
[40] P. Thalau, T. Ritz, K. Stapput, R. Wiltschko, and W. Wiltschko, Naturwissenschaften 92, 
86 (2005). 

[41] N. Keary, T. Ruploh, J. Voss, P. Thalau, R. Wiltschko, W. Wiltschko, and H.-J. Bischof, 

Frontiers in Zoology 6 (2009). 
[42] T. Ritz, P. Thalau, J. B. Phillips, R. Wiltschko, and W. Wiltschko, Nature 429, 177 (2004). 
[43] T. Ritz, R. Wiltschko, P. Hore, C. T. Rodgers, K. Stapput, P. Thalau, C. R. Timmel, and 

W. Wiltschko, Biophys. J. 96, 3451 (2009). 
[44] R. Wiltschko, K. Stapput, T. Ritz, P. Thalau, and W. Wiltschko, HFSP J. 1, 41 (2007). 
[45] I. A. Solov'yov, D. E. Chandler, and K. Schulten, Plant Signaling & Behavior 3, 676 (2008). 
[46] F. Cintolesi, T. Ritz, C. Kay, C. Timmel, and P. Hore, Chem. Phys. 294, 707 (2003). 
[47] K. Maeda, K. B. Henbest, F. Cintolesi, I. Kuprov, C. T. Rodgers, P. A. Liddell, D. Gust, 

C. R. Timmel, and P. J. Hore, Nature 453, 387 (2008). 
[48] M. Ahmad and A. R. Cashmore, Plant Mol. Biol. 30, 851 (1996). 
[49] A. R. Cashmore, J. A. Jarillo, Y.-J. Wu, and D. Liu, Science 284, 760 (1999). 
[50] A. Sancar, Chem. Rev. 103, 2203 (2003). 
[51] M. Ahmad and A. Cashmore, Nature 366, 162 (1993). 
[52] C. Lin and T. Todo, Gen. Biol. 6, 220.1 (2005). 
[53] C. L. Partch and A. Sancar, Meth. Enzym. 393, 726 (2005). 
[54] J. M. Christie and W. R. Briggs, J. Biol. Chem. 276, 11457 (2001). 

[55] M. Liedvogel, K. Maeda, K. Henbest, E. Schleicher, T. Simon, C. R. Timmel, P. Hore, and 

H. Mouritsen, PLoS ONE 2, ell06 (2007). 
[56] T. Biskup, E. Schleicher, A. Okafuji, G. Link, K. Hitomi, E. D. Getzoff, and S. Weber, 

Angew. Chem. Int. Ed. Engl. 48, 404 (2009). 



32 



[57] Y.-T. Kao, C. Tan, S.-H. Song, N. Oztiirk, J. Li, L. Wang, A. Sancar, and D. Zhong, J. Am. 

Chem. Soc. 130, 7695 (2008). 
[58] T. R. Prytkova, D. N. Beratan, and S. S. Skourtis, Proc. Natl. Acad. Sci. USA 104, 802 

(2007). 

[59] K. B. Henbest, K. Maeda, P. J. Hore, M. Joshi, A. Bacher, R. Bittl, S. Weber, C. R. Timmel, 

and E. Schleicher, Proc. Natl. Acad. Sci. USA 105, 14395 (2008). 
[60] M. Byrdin, A. P. Eker, M. H. Vos, and K. Brettel, Proc. Natl. Acad. Sci. USA 100, 8676 

(2003) . 

[61] Y. M. Gindt, E. Vollenbroek, K. Westphal, H. Sackett, A. Sancar, and G. T. Babcock, 

Biochem. 38, 3857 (1999). 
[62] C. Aubert, P. Mathis, A. P. Eker, and K. Brettel, Proc. Natl. Acad. Sci. USA 96, 5423 

(1999). 

[63] C. Aubert, M. H. Vos, P. Mathis, A. P. Eker, and K. Brettel, Nature 405, 586 (2000). 

[64] H. Mouritsen, U. Janssen-Bienhold, M. Liedvogel, G. Feenders, J. Stalleicken, P. Dirks, and 

R. Weiler, Proc. Natl. Acad. Sci. USA 101, 14294 (2004). 
[65] A. Moller, S. Sagasser, W. Wiltschko, and B. Schierwater, Naturwissenschaften 91, 585 

(2004) . 

[66] M. Liedvogel and H. Mouritsen, J. R. Soc. Interface 7, S147 (2010). 

[67] C. NieBner, S. Denzau, J. C. Gross, L. Peichl, H.-J. Bischof, G. Fleissner, W. Wiltschko, and 

R. Wihschko, PLoS ONE 6, e20091 (2011). 
[68] H. Mouritsen, G. Feenders, M. Liedvogel, K. Wada, and E. Jarvis, Proc. Natl. Acad. Sci. 

USA 102, 8339 (2005). 

[69] M. Liedvogel, G. Feenders, K. Wada, N. Troje, E. Jarvis, and H. Mouritsen, Eur. J. Neurosci. 
25, 1166 (2007). 

[70] G. Feenders, M. Liedvogel, M. Rivas, M. Zapka, H. Horita, E. Kara, K. Wada, H. Mouritsen, 

and E. Jarvis, PLoS ONE 3, el768 (2008). 
[71] D. Heyers, M. Manns, H. Luksch, O. Giintiirkiin, and H. Mouritsen, PLoS ONE 2, e937 

(2007). 

[72] C. M. Hein, M. Zapka, D. Heyers, S. Kutzschbauch, N.-L. Schneider, and H. Mouritsen, J. 

R. Soc. Literface 7, S227 (2010). 
[73] W. Wiltschko, J. Traudt, O. Giintiirkiin, H. Prior, and R. Wiltschko, Nature 419, 467 (2002). 



33 



[74] C. M. Hein, S. Engels, D. Kishkinev, and H. Mouritsen, Nature 471, Ell (2011). 
[75] W. Wiltschko, J. Traudt, O. Giintiirkiin, H. Prior, and R. Wiltschko, Nature 471, E12 
(2011). 

[76] K. Stapput, O. Guntiirkiin, K.-P. Hoffmann, R. Wiltschko, and W. Wiltschko, Cur. Biol. 20, 
1259 (2010). 

[77] J. -P. Bouly, E. Schleicher, M. Dionisio-Sese, F. Vandenbussche, D. Van Der Straeten, 
N. Bakrim, S. Meier, A. Batschauer, P. Galland, R. Bittl, et al., J. Biol. Chem. 282, 9383 
(2007). 

[78] B. Giovani, M. Byrdin, M. Ahmad, and K. Brettel, Nat. Struct. Biol. 10, 489 (2003). 

[79] N. Hoang, E. Schleicher, S. Kacprzak, J.-P. Bouly, M. Picot, W. Wu, A. Berndt, E. Wolf, 

R. Bittl, and M. Ahmad, PLoS Biol. 6, 1559 (2008). 
[80] J. Shirdel, P. Zirak, A. Penzkofer, H. Breitkreuz, and E. Wolf, Chem. Phys. 352, 35 (2008). 
[81] C. A. Brautigam, B. S. Smith, Z. Ma, M. Palnitkar, D. R. Tomchick, M. Machius, and 

J. Deisenhofer, Proc. Natl. Acad. Sci. USA 101, 12142 (2004). 
[82] R. Pokorny, T. Klar, U. Hennecke, T. Carell, A. Batschauer, and L.-O. Essen, Proc. Natl. 

Acad. Sci. USA 105, 21023 (2008). 
[83] R. Banerjee, E. Schleicher, S. Meier, R. M. Viana, R. Pokorny, M. Ahmad, R. Bittl, and 

A. Batschauer, J. Biol. Chem. 282, 14916 (2007). 
[84] T. Langenbacher, D. Immeln, B. Dick, and T. Kottke, J. Am. Chem. Soc. 131, 14274 (2009). 
[85] A. Zeugner, M. Byrdin, J.-P. Bouly, N. Bakrim, B. Giovani, K. Brettel, and M. Ahmad, 

J. Biol. Chem. 280, 19437 (2005). 
[86] A. Berndt, T. Kottke, H. Breitkreuz, R. Dvorsky, S. Hennig, M. Alexander, and E. Wolf, 

J. Biol. Chem. 282, 13011 (2007). 
[87] N. Oztiirk, S. Song, C. Selby, and A. Sancar, J. Biol. Chem. 283, 3256 (2008). 
[88] J. B. Phillips, P. E. Jorge, and R. Muheim, J. R. Soc. Interface 7, S241 (2010). 
[89] B. Liu, H. Liu, D. Zhong, and C. Lin, Cur. Opin. Plant Biol. 13, 578 (2010). 
[90] J. Brazard, A. Usman, F. Lacombat, C. Ley, M. M. Martin, P. Plaza, L. Mony, M. Heijde, 

G. Zabulon, and C. Bowler, J. Am. Chem. Soc. 132, 4935 (2010). 
[91] L Chaves, R. Pokorny, M. Byrdin, N. Hoang, T. Ritz, K. Rettel, L.-O. Essen, G. T. van der 

Horst, L. Batschauer, and M. Ahmad, Annu. Rev. Plant Biol. 62, 335 (2011). 
[92] K. E. O'Day, PLoS Biol. 6, 1359 (2008). 



34 



[93] M. Byrdin, V. Sartor, A. P. Eker, M. H. Vos, C. Aubert, K. Brettel, and P. Mathis, Biochim. 
Biophys. Acta 1655, 64 (2004). 



[94 
[95 
[96 

[97; 

[98 
[99 
[100 
[101 
[102 
[103 
[104 
[105 

[106 
[107; 
[108 
[109 

[110 
[111 

[112 
[113 
[114 



T. Kottke, A. Batschauer, M. Ahmad, and J. Heberle, Biochem. 45, 2472 (2006). 
K. Schulten, J. Chem. Phys. 82, 1312 (1985). 

H. J. Hogben, O. Efimova, N. Wagner-Rundell, C. R. Timmel, and P. More, Chem. Phys. 
Lett. 490, 118 (2009). 

W. Schmidt, Optical Spectroscopy in Chemistry and Life Sciences, an Introduction (Wiley- 
VCH, 2005). 

R. J. Gegear, A. Cassehnan, S. Waddell, and S. M. Reppert, Nature 454, 1014 (2008). 

L. E. Foley, R. J. Gegear, and S. M. Reppert, Nature (2011). 

O. Efimova and P. More, Biophys. J. 94, 1565 (2008). 

P. Semm and R. Reason, Brain Res. Bull. 25, 735 (1990). 

T. Alerstam, J. Exp. Biol. 130, 63 (1987). 

A. Lednor and C. Walcott, Beh. Ecol. Sociobiol. 130, 3 (1988). 

U. Munro, J. Munro, and J. Phillips, Naturwissenschaften 84, 26 (1997). 

D. Heyers, M. Zapka, M. HofFmeister, J. Wild, and H. Mouritsen, Proc. Natl. Acad. Sci. USA 

107, 9394 (2010). 

D. Kishkinev, N. Chernetsov, and H. Mouritsen, Auk 127, 773 (2010). 

W. Wiltschko and R. Wiltschko, J. Comp. Physiol. A 184, 295 (1999). 

W. Wiltschko, R. Wiltschko, and U. Munro, Naturwissenschaften 87, 366 (2000). 

W. Wiltschko, U. Munro, H. Ford, and R. Wiltschko, Proc. R. Soc. Lond. B 270, 2133 

(2003). 

J. R. Lakowicz, Limological Analysis (Springer, New York, 2006). 

A. Bartels, A. Hudert, C. Janke, T. Dekorsy, and K. Kohler, Appl. Phys. Lett. 88, 04117 
(2006). 

R. Berera, R. van Grondelle, and J. T. M. Kennis, Photosynth. Res. 101, 105 (2009). 
C. Lin and D. Shalitin, Annu. Rev. Plant Biol. 54, 469 (2003). 

T. Ritz, T. Yoshii, C. Helfrich-Foerster, and M. Ahmad, Commun. &: Integr. Biol. 3, 21 
(2010). 



